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1.0  INTRODUCTION 


The  problem  of  damage  assessment  of  interface  circuits  by  EMP- 
generated  pulses  is  of  much  interest  to  the  E-3A  Program,  as  well  as 
other  aircraft  systems  which  must  satisfy  EMP  nuclear  requirements. 
From  network  analysis  of  the  interface  circuits,  damage  assessment  is 
accomplished  by  evaluating  whether  the  voltage  pulses  which  appear 
across  the  susceptible  components  will  damage  the  latter.  References 
1 and  2 collectively  provide  a good  discussion  of  the  damage  mecha-r 
nisms  for  various  components. 

The  purpose  of  this  investigation  is  to  discuss  some  analytical 
techniques  which  can  be  used  for  evaluating  permanent  damage  of  semi- 
conductor devices  which  are  contained  in  interface  circuits.  The 
results  are  applicable  to  the  study  of  failure  of  semiconductor  June- 

(3) 

tions,  as  originally  proposed  by  Wunsch  . In  his  model,  the 
destruct  mechanism  is  assumed  to  result  from  changes  in  the  junction 
parameters  due  to  the  high  temperatures  produced  locally  within  the 
junction  area.  These  hot  spots  can  be  produced  for  both  forward  and 
reverse  voltage  conditions.  Thus,  the  evaluation  of  component  suscep- 
tibility reduces  to  determining  whether  the  temperature  increase  due 
to  power  dissipation  is  sufficient  to  cause  damage. 

Damage  assessment,  therefore,  depends  upon  the  power  delivered 
to  a component,  which  in  turn  depends  on  the  signal  delivered  to  that 
component  through  a network  of  circuitry.  The  complexity  of  the  cal- 
culation can  vary  considerably.  Some  cases  are  simple  enough  for  hand 

(4) 

analysis  while  others  require  computer  codes  . Because  computer 
analysis  costs  can  be  very  high  in  a system  containing  a large  number 
of  components,  the  preliminary  hand  analysis/screening  of  interface 
circuits  is  an  important  step  in  the  hardness  evaluation  process.  An 
analytical  assessment  not  only  provides  a reasonable  assessment  of 
circuit  hardness,  but  also  identifies  those  circuits  requiring  a more 
vigorous  assessment* 
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In  this  study,  we  provide  a discussion  of  some  of  the  physical 
ideas  and  analytical  concepts  which  are  involved  in  hand  analysis  of 
component  damage.  It  is  not  unlikely  that  a portion  of  what  is  pre- 
sented here  exists  in  other  reports,  but  this  information  does  not 
seem  to  have  been  conveniently  coordinated.  Therefore,  this  document 
also  provides  the  reader  with  a general  background  of  the  problem. 

The  remainder  of  this  report  is  organized  as  follows:  Section 

2.0  deals  with  the  extension  of  Wunsch’s  analysis  to  finite  size  p-n 
junctions.  The  subsequent  section  discusses  various  aspects  pertain- 
ing to  the  analytical  behavior  of  junction  temperature  including  the 
effect  of  interpulse  temperature  cooling;  and  the  response  of  inter- 
face circuits  to  damped  sine  inputs.  Concluding  remarks  are  rendered 
in  Section  4.0. 
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2.0  DERIVATION  OF  MAXIMUM  TEMPERATURE  IN  P-N  JUNCTION  DIODE 


in  a p-n  junction  in  an  infinite  medium.  The  purpose  of  this  section 
of  the  report  is  to  extend  the  analysis  to  p-n  junctions  of  finite 
size.  Using  the  developed  equations,  one  can  readily  assess  the  pulse 
length  conditions  for  which  the  finite  size  of  the  junction  should  be 
taken  into  account  even  in  a hand  screening  analysis,* 

Figure  1 shows  a pictorial  representation  of  a p-n  junction  diode. 
For  analysis  purposes,  the  p and  n sides  are  taken  to  be  of  equal 
dimension,  L/2,  with  the  junction  located  at  L/2, 

When  a voltage  pulse  is  applied  to  a semiconductor  device  in  the 
reverse  direction,  the  principal  voltage  drop  is  across  the  junction. 

If  we  assume  that  the  junction  is  infinitesimally  thin,  then  the 
spatial  distribution  of  volumetric  heating  (watts/cm3)  is  given  by 


where  P(t)  is  the  time-dependent  total  power  delivered  to  the  device, 
A is  the  cross  sectional  area  of  the  junction  and  6(x)  is  the  Dirac 
delta  function.  The  general  one-dimensional  heat  diffusion  equation 


(2.1) 


is  given  by^: 


^(PCT)  ■ £ <k  g ) ♦ «*.t) 


(2.2) 


where: 


p ■ density  ( ^ ) 
cm  ^ 


C = specific  heat 


^ 1 oules  j 
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T = temperature  (°K) 

k = thermal  conductivity  ( ) 

cm-°K 

x = linear  dimension  (cm) 


Junction  (area=A) 


L 

2 


L 

2 


x=0 


x=L 


Figure  1.  Diagram  of  p-n  Junction  Diode 
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Wunsch  points  out  that  k is  a sensitive  function  of  temperature, 
e.g.,  ranging  from  1.56  watts/cm-°K  at  300°K  to  0.310  at  1000°K,  so 
that  strictly  speaking  its  explicit  dependence  on  T should  be  taken 
into  account  in  Equation  (2,2).  However,  for  mathematical  conveni- 
ence, we  shall  approximate  k by  a suitable  time-weighted  average  over 

(3) 

the  temperature  range  of  interest 

Equation  (2.2)  will  be  solved  subject  to  the  boundary  conditions 


T(x*=0)  - T 

a 

T(x=L)  = T 

a 


(2.3) 


where  T is  the  ambient  temperature  which  may  also  be  taken  as  the 
a 

initial  temperature.  Without  loss  of  generality,  we  arbitrarily  mea- 
sure the  temperature  with  respect  to  T . We,  thus,  make  the  substitu- 

a 

tion  T-+-T-T  so  that  the  new  boundary  conditions  become  T(x=0)  = 

cl 

T(x*L)=0.  Using  the  aforementioned  change  of  temperature  reference, 
the  solution  of  Equation  (2.2)  is  obtained  by  using  an  orthogonal 
series  expansion  for  T(x,t).  We  let 


T(x,t) 


00 


(t)  sin 


mrx 


(2.4) 


q(x,t) 


00 


(t)  sin 


mrx 


(2.5) 


where  the  orthogonal  functions  sin 


nTTX 


satisfy  the  condition: 


/,  mrx  mux  L / 

sin  — — sin  — — dx  ® — / i 

L L it  J 


sin  nu  sin  mu  du 


(2.6) 


if  m = n 


0 if  m ^ n 
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The  functions  Q (t)  are  determined  from  the  equation: 
n 


Li 

/ 


sin  q (x,t)  dx  = — Q (t) 

L z n 


(2.7) 


which  using  Equation  (2.1)  yields; 


Q„(t) 


2 

L 


A 


sin 


nTr 


(2-8) 


For  mathematical  convenience,  let  us  call 


. nn 
sin  — 


n 


0 for  n even 


(2-9) 


±1  for  n odd 


Substituting  Equations  (2,4)  and  (2.5)  into  Equation  (2.1)  and  equating 
equal  coefficients  of  sin  gives: 


PC 


dB 
n 

dt 


-k  (S£)  B + Q (t) 
L n n 


(2.10) 


Since  C = 0 for  n even.  Equation  (2.10)  has  meaning  for  odd  n,  a 
n 

fact  we  will  later  use  in  summing  up  the  infinite  series.  Dividing 
through  by  pC  gives  the  following  simplified  version  of  Equation 
(2.10); 


dB 
n 

dt 


-X 


n 


B + W (t) 
n n 


(2.11) 
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where 


(2.12) 

(2.13) 

The  initial  conditions  for  all  B are: 

n 

B (t=0)  = 0 , (2.14) 

n 


k<f>‘ 


w 


pc 


2P(,t) 

pCAL 


since  it  is  assumed  initially  that  the  junction  is  at  uniform  tempera- 
ture T . The  solution  to  Equation  (2.11),  subject  to  the  condition 
a 

of  Equation  (2.14)  is  given  by: 

-Ant  r Ant 1 

B (t)  = e / W (t1)  e dt1  (2.15) 

n J n 

0 


On  physical  grounds,  it  is  apparent  that  the  maximum  temperature 
will  occur  at  the  junction,  where  x = Let  us  define  this  junction 

temperature  as  T.  We  have: 

OO  00 

T(t)  - T(x=|,t)  = X)  Bn(t)  Sin  It  2 = HBn(t)  cn  (2,16) 

n=l  n=l 

(n  odd)  (n  odd) 


where  C is  given  by  Equation  (2.9).  Substituting  Equation (2. 15) 
n 


into  Equation  (2.16)  gives 


pCAL 


/P<t,,(£ 

0 ' n=l 


-Xn(t-t') 


dt1 


(2.17) 


1 n=l 
(n  odd) 
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Since  according  to  Equation  (2*9),  C = ±1,  it  then  follows 

2 n 
that  = 1 so  that  the  infinite  summation  in  Equation  (2,17) 

becomes : 


S 


a y.  "Xn*t"t  p 2 
' L.  e Cn 

n=l 

(n  odd) 


-Xj(t-t f)n2 


n=l 

(n  odd) 


oo 


n=l 

(n  odd) 


(2.18) 


where 


(2.19) 


a = (t-t') 


(2.20) 


S is  computed  in  the  following  way.  We  first  convert  Equation 
(2.18)  to  an  unrestricted  series  through  the  substitution: 


2m  + 1 


(2.21) 


where  the  index  m ranges  from  m=0  to  m»  00 . S then  becomes: 

OO  OO 

-a(2nrt-l)2  _ V"  -4a(nH-l/2)2 


£ * 

m=0 


Yj  6 

m=0 


(2.22) 


If  a is  always  very  small,  as  for  example  would  be  the  case  in 
an  infinite  medium  for  which  L-*»  (cf.  Equation  2.19),  then  it  is  clear 
that  each  succeeding  term  in  the  summation  will  only  differ  a small 
amount  from  its  predecessor.  In  this  case,  the  discrete  summation  can 
readily  be  approximated  by  /dm  . We  have 

00 

-4ay2 


Lim  S 
a-K) 


■/ 

0 


c~4a  (nri-1/2)  2dm 


■/ 

1/2 


dy 


(2.23) 
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Making  the  substitution 


w2  = 4**/ 


gives 


Lira  S 
a-K) 


Since 


we  have 

/rr  >/n  L fic  /pc  L 

Lim  S = * - ■ = — 

a-*-0  4 \Ta  4 \fk  tt  /t-t * 4 \/kTr(t-t*) 


(2,24) 


(2.25) 


(2.26) 


(2.27) 


Substituting  Equation  (2.27)  into  Equation  (2,17)  then  gives: 


T(t) 


pCAL 


/ PCL  /* 

4 ^ 


P(t')  dt' 

/^tT 


1/2  t. 

2 WkpC  / J 


dt' 

1/2 


(2.28) 


By  making  the  appropriate  notational  substitutions,  and  generalizing 
the  results  of  reference  3 to  the  case  of  an  arbitrary  time  history  of 
P(t),  the  reader  can  convince  himself  that  Equation  (2.28)  is  exactly 
equal  to  the  result  derived  by  Wunsch.  For  example,  if  we  assume: 

P(t)  - P : 0 < t < t (2.29) 

m o — — s 


then  Equation  (2.28)  yields  the  well-known  formula: 
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1/2 


T(t) 

*°\< 

H 

( 1 ) 

v irkpC  / 

1/2 

t1/2  : 0 < t < t 

— — s 

(a) 

A 

T 

max 

P 

- f 

( 1 ) 

\ irkpC  ' 

1/2 

i t 1/2 

s 

(b) 

(2.30) 

Examination  of  Equations  (2.19)  and  (2.20)  show  that  in  general 
a will  be  small  whenever 

« 1 (2.31) 

or  equivalently 

t«td-f(|)  (2.32) 

where  t,  is  recognized  as  the  diffusion  time.  If  the  pulse  duration 
d 

is  t then  Equation  (2.32)  shows  that  for  t <<  t,  the  result  given 
P P d 

by  Equation  (2.28)  provides  an  accurate  estimate  of  the  junction  tem- 
perature. 

On  the  other  hand,  when  a becomes  comparable  to  or  greater  than 
unity,  the  approximation  of  the  infinite  series  by  an  integral  becomes 
questionable  and  an  alternate  method  must  be  found  for  evaluating 

A 

T(t).  Fortunately,  this  can  be  done  by  performing  a term-by-term  time 
integration  in  Equation  (2.17).  Making  the  substitution 

2 

X = 1 n (2.33) 

n ± 

and  using  Equation  (2.29)  for  P(t)  in  Equation  (2.17)  for  illustra- 
tive purposes,  we  obtain: 


V. 


10 


oo 


n=l 

(n  odd) 


(2.34) 


OO 


(n  odd) 


Theoretically  speaking,  Equation  (2.34)  is  correct  for  all  times, 
but  the  infinite  series  has  poor  convergence  properties  in  the  limit 
of  small  \±t.  On  the  other  hand,  the  series  will  converge  rapidly 
when  A^t  1 which  corresponds  to  the  regime  where  t _>  t, . In  this 
time  regime,  the  exponential  terms  will  decay  rapidly  with  T(t)  reach- 
ing an  asymptotic  value  given  by 


T (t)  could  have  also  been  derived  from  a steady-state  solution  of 
max 

Equation  (2.2). 


n=l  11 
(n  odd) 


(2.35) 


From  reference  6,  we  can  deduce 


OO 


oo 


(2.36) 


(n  odd) 


which  when  substituted  into  Equation  (2,35)  yields 


P L 


o 


(2.37) 


4AK 
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It  is  worthwhile  to  mention  that,  for  long  pulse  lengths,  as 
defined  by  the  condition  t >>  t^,  a steady-state  solution  is  meaning- 
ful. In  this  case,  the  explicit  temperature  dependence  of  the  conduc- 
tivity can  be  taken  into  account  in  the  computation,  since  in  the  p 
and  n regions  we  deal  with  the  equation 


0 


_3_ 

3x 


(2.38) 


This  equation  integrates  to: 

T 

k (T')  dT* 


/' 


Cx 


(2,39) 


where  C is  a constant.  By  straightforward  application  of  the  boundary 
conditions,  the  spatial  distribution  of  temperature  can  be  determined. 

A 

In  summary.  Equation  (2.34)  provides  a means  of  calculating  T(t) 
for  all  times  and  finite  size  p-n  junctions.  For  the  limiting  case  in 
which  (t/tjj)  <<  1,  Wunsch's  original  result  is  valid.  At  the  other 

A 

extreme,  when  (t/t^)  >>  1,  the  asymptotic  bound  on  T(t)  should  be  used. 
For  intermediate  times,  a finite  number  of  terms  from  Equation  (2.34) 
would  probably  provide  an  adequate  approximation. 
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3.0  ANALYTIC  BEHAVIOR  OF  JUNCTION  TEMPERATURE 


The  purpose  of  this  section  is  to  discuss  various  analytical  pro- 
perties of  the  time  dependence  of  the  junction  temperature  as  a func- 
tion of  power  pulse  shape.  The  results  may  be  useful  in  a screening 
assessment  of  damage  susceptibility. 


3.1  Alternate  Form  of  Damage  Criterion 


For  the  purposes  of  this  discussion,  we  shall  assume  that  the 

junction  temperature  is  given  by  Equation  (2.28).  For  a square  wave 

pulse  of  power  Pq , and  duration  tg,  the  maximum  temperature  is 

given  by  Equation  (2.30b).  If  Tm  is  the  temperature  at  which  damage 

occurs,  the  following  well-known  relationship  holds  between  T * P , 

m o 

and  t : 
s 


(3.1) 


where  K is  the  damage  constant  defined  by  the  equation 


K = A /irkpC  Tm  (3.2) 

For  this  special  case  of  a square  pulse,  the  total  energy  delivered  to 
the  device  is: 


E - P t 
s os 

which  gives  the  equivalent  relationship 


E 

s 


Kt 


1/2 


s 


(3.3) 


(3.4) 
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Using  Equation  (3.2),  we  can  recast  Equation  (2,28)  in  the  mathe- 
matically convenient  form 


T(t) 


1 GM 


(3.5) 


where 


G(t) 


P(t;) 


(t-t') 


1/2 


dt' 


(3.6) 


Equations  (3.5)  and  (3.6)  provide  an  alternate  way  of  evaluating 

damage.  Since  damage  will  occur  when  T(t)  exceeds  T , the  condi- 

m 

tion  for  damage  to  occur  is  obtained  from  the  inequality 


G(t)  >_  2K  (3.7) 

The  foregoing  expression  for  damage  appears  to  be  a new  way  of 

relating  pulse  shape  to  component  susceptibility.  If  is  the 

maximum  value  of  G(t)  in  the  interval  0 < t < t , where  t is  the 

~~  ~ P P 

pulse  duration,  then  damage  will  occur  when 
G 

max 

Even  for  pulses  of  finite  duration,  it  is  clear  that  the  maximum 
of  G(t)  need  not  occur  at  the  end  of  the  pulse  interval,  so  that 
damage  formulas  which  are  based  on  total  energy  deposition 

t 

p(t')  dt’ 

o 

may  lead  to  uncertainties  in  damage  assessment.  However,  these  uncer- 
tainties will  not  be  very  large  for  the  commonly-assumed  single  cycle 


max 


!/ 


PCt*) 


(t-t') 


1/2 


dt' 


2K 


(3.8) 
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waveforms  (e.g.,  one  cycle  of  a sine  wave).  For  example,  for  a limited 
number  of  cases  evaluated,  it  was  found  that  about  75%  of  the  pulse 
energy  contributed  to  raising  the  junction  temperature  to  its  maximum 
value . 

Equation  (3.6)  may  be  useful  for  determining  the  maximum  tempera- 
ture for  complex  waveforms,  and  also  for  more  accurately  assessing 
(cf.  Section  3.3)  the  damped  sine  case  than  in  reference  2.  In  addi- 
tion, it  provides  an  easy  means  of  computing  the  interpulse  cooling 
rate  for  analyses  of  repeated  pulses  (cf.  Section  3.2), 

3.2  Interpulse  Temperature  Cooling 

The  purpose  of  this  section  is  to  assess  the  temperature  decrease 
between  pulses.  For  mathematical  simplicity,  we  consider  the  pulse 
train  shown  in  Figure  2.  Using  the  structure  for  P(t)  shown  in  this 
figure. 


P 


P(t) 


o 


t > 


Figure  2.  Power  Pulse  Train 


the  temperature  is  determined  from  the  equation: 


0 


(3.9) 


In  the  range  0 <.  t _<  T,  the  function  G(t)  becomes: 


1/2 

0 < t < T:  G(t)  = Gn(t)  = 2P  tA/ 

— — j.  o 


(3.10) 
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while  in  the  range  T ^ t £ 2T,  we  have: 

T < t < 2T:  G(t)  = G_(t)  = 2P  [t1^2-  (t-T)1^2]  (3.11) 

— z o 


The  function  G2^  gives  the  temperature  response  following  a square 
wave  power  impulse  function  and  is  actually  valid  for  all  times  greater 
than  T.  For  t >>  T,  the  function  G2^  can  aVVrox±mated  by  the 
expansion: 


G2  (t) 


2P  [t1/2  - t1/2U  - V/2] 
o t 


2P  [t1'2  - t1/2(l 
o 


1 1 - 1 (I))?  J 
2 t 8 V'  1 


(3.12) 


T 1 T ' 
(P  -)  + f P (-) 
v o t A o t 


The  asymptotic  expansion  can,  of  course,  be  directly  obtained  from 
Equation  (3.9)  by  neglecting  t'  with  respect  to  t in  the  denomina- 
tor of  the  integral.  For  times  t long  compared  to  the  pulse  width, 

A 

and  for  any  shape  pulse,  T(t)  becomes: 


T(t) 


T 

m 

2K 


P(t’)  dt’ 


T E 
m _g_ 


2K  tl/2 


(3.13) 


where 


E 

P 


t 


P(t’)  dt' 


energy  delivered  in  pulse. 


(3.14) 


Figure  3 shows  a plot  of  G(t)  (normalized  to  2P  T = 1)  for  the 

o 

first  pulse  of  the  pulse  train;  the  function  G£(t)  is  shown  for  all 
time.  An  interesting  feature  of  G2(t)  is  its  relatively  slow  decay 
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Figure  3.  Plot  of  G(t)  for  First  Pulse 


Figure  4.  Plot  of  G(t)  for  Three  Pulses 
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with  time.  This  is  a reflection  of  Wunsch's  physical  model  which 
limits  the  amount  of  cooling  for  an  infinite  system. 

The  function  G(t)  at  times  t = (2n+l)T  is  given  by 

n 

G[(2n+1)T]  = G1(T)  + ^ G2  [(2n’+l)T]  (3.15) 

n ,=1 

Figure  4 shows  a plot  of  G(t)  for  three  pulses.  For  large  n,  the 
function  G2(t)  can  be  approximated  by  the  first  term  in  Equation 
(3.12)  with  the  normalization  PqT  = 1/2.  Thus,  for  large  n we  have: 

G2t(2n+1)T)  = J (2n+l)T  = 2(2n+l)  (3.16) 

which  shows  that  the  series  of  Equation  (3.15)  is  logarithmically  diver- 
gent. The  upshot  of  this  discussion  is  that  it  is  precarious  to  make 
computations  of  damage  assessment  based  upon  one  cycle  of  a series  of 
pulses.  For  the  Wunsch  model,  temperature  quenching  appears  to  be 
generally  insufficient. 

The  foregoing  discussion  is  useful  for  analysis  of  exponentially 

damped  periodic  pulses  (which  includes  the  damped  sine  case) . Let 

T be  the  period  of  the  oscillation,  and  Y the  damping  constant, 
c 

The  power  waveform  is  expressed  as 

P(t)  = Ae"Yt  W(t)  (3.17) 

where  W(t)  is  defined  as  the  periodic  part  which  satisfies  the 
condition 


W(t+T  ) « W(t) 

c 


(3.18) 


In  order  for  Equation  (3.17)  to  have  any  meaning,  it  is  clear  that 
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YT  <<  1.  If  damping  is  now  applied  to  the  pulse  waveform  in  Figure  2 

(T  = 2T) , with  yt  being  considered  to  be  relatively  constant  between 

jaT  and  (nrKL)T  , and  approximated  by  the  value  mT  during  the  inter- 
c c 

val,  the  the  corresponding  behavior  of  G( t)  at  the  odd  intervals  of 

T would  be: 

n 

G[(2n+1)T]  = e"2nYTG1(T)  + ^ e 2 (n_n' ^YTG2 [ (2n'+l)T]  (3.19) 

n'=l 


For  this  situation,  it  is  not  clear  during  what  cycle  the  maximum 
temperature  will  occur.  Equation  (3.19)  does,  however,  provide  a means 
of  readily  evaluating  T[(2n+1)T],  and  thereby  determining  the  number  of 
cycles  which  must  be  incorporated  into  the  calculation. 

3. 3 Response  of  Interface  Circuits  to  Damped  Sine  Inputs 

It  has  been  found  experimentally  and  shown  theoretically  that  the 
voltage  input  to  interface  circuits  is  frequently  of  the  damped  sine 
type.  Such  a voltage  pulse  is  described  by  the  form: 

V(t)  = VQ  e“Ytsin  0)t  (3.20) 

where  in  order  for  Equation  (3.20)  to  have  meaning  we  must  have 

^ » 1 (3.21) 

Voltages  of  the  type  given  by  Equation  (3.20)  may  give  rise  to  various 
forms  of  power  dissipation  in  devices  depending  upon  whether  failure 
takes  place  in  the  forward  or  backward  direction.  Some  of  the  typical 
power  waveforms  are  discussed  in  references  2 and  4.  Figure  5 shows 
some  representative  power  waveforms  resulting  from  damped  sine  inputs. 
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If  the  interface  circuit  is  purely  resistive,  then  the  power  dis- 
sipation in  a particular  device  will  be  of  the  periodic  type,  not 

unlike  those  shown  in  Figure  5.  On  the  other  hand,  if  failure  takes 

place  in  the  reverse  direction,  the  length  of  time  for  conduction  will 

(2) 

be  voltage  dependent  and  the  power  dissipation  will  not  be  strictly 
periodic.  For  a damped  sine  input,  we  will  eventually  reach  the  con- 
dition where  the  reverse  voltage  falls  below  the  breakdown  voltage  and 
conduction  in  the  negative  direction  ceases. 

When  inductive  and  capacitive  effects  have  to  be  taken  into 
account  in  the  interface  circuit  model,  it  is  not  clear  that  power 

dissipation  will  be  cyclic  since  the  basic  problem  is  a transient  one. 

Circuit  conditions  in  the  second  and  subsequent  cycles  of  the  input 
voltage  may  be  different  from  the  first,  depending  on  the  values  of 
the  reactive  elements.  Therefore,  in  some  cases,  a time-dependent 
calculation  would  appear  to  be  the  only  sensible  method  of  solution. 

However,  there  do  exist  several  instances  where  effects  arising 
from  non-zero  initial  conditions  quickly  dampen  out,  thus  permitting 
analyses  of  each  cycle  separately.  The  following  comments  provide  a 
discussion  of  damage  assessment  under  these  conditions. 

The  first  logical  step  in  damage  assessment  is  to  approximate 
V(t)  by 


V(t)  = V sin  ut  (3.22) 

o 

for  the  first  few  cycles  subject  to  the  condition  yt  < 1.  In  view 
of  the  fact  that  (w/y)  » 1,  this  time  regime  may  include  many  cycles, 
If  damage  is  found  to  occur  in  the  first  few  cycles,  it  is  then 
unnecessary  to  include  the  damping.  In  performing  such  an  analysis, 
it  would  of  course  be  necessary  to  use  the  results  of  the  previous 
section  to  account  for  temperature  decay. 
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i 


P(t) 


Figure  5.  Representative  Power  Waveforms 
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However,  additional  insight  into  the  problem  can  be  found  with 
use  of  the  following  observation:  In  the  periodic  response  range 

(characterized  by  Equation  3.22),  the  power  P(t)  across  any  device 
can  be  cast  in  the  form: 

00 

P(t)  - P + y P sin  mot  (3.23) 

ave  - n 

n=l 

where  P is  the  average  of  P(t)  and  P are  the  Fourier  coeffi~ 
ave  n 

cients.  Equation  (3.23)  would  be  valid  in  the  regime  yt  < 1 which, 
as  previously  indicated,  could  include  the  range  cot  > 1,  If  Equation 
(3.23)  is  inserted  in  the  expression  for  G(t),  we  obtain: 


^ = 2p  t1/2  + y> 

1/2  ave  L-j  n 

n=l 


t 

/ 


sin  ncot' 


(t-t’) 


1/2 


dt'  (3.24) 


G(t) 


„ f lit 

l (t-f) 


Working  out  the  details  for  the  integration  in  Equation  (3.24),  we 
obtain: 


t 

i (« . f j2 


dt' 


where 


(4"  r 


sin  ncot  f-  (t) 
In 


cos  ncot  f0  (t) 
zn 


(3.25) 


f2n(t) 


ncot 

■/ 


cos  1 

1/2 


dA 


(3.26) 


ncot 

/ 


sin  A 

,1/2 


dA 


As  per  the  previous  discussion,  the  foregoing  set  of  equations 
are  of  interest  in  the  time  regime  encompassing  many  cycles,  but  not 
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long  enough  to  include  damping.  In  the  regime  where  cot  >>  1 (e.g., 

cot=13  for  2 cycles),  it  can  be  shown^  that 


ncot 

/ 


sin  X 

,1/2 


dX 


mot 

I 


cos  X 

1/2 


dX 


■N/f 


(3.27) 


in  which  case  I (t)  becomes: 
n 

1/2 


Xn(t)  = yp  (i)  tSin  nt0t  " COS  nU)t]  ^8in  (mot  ‘ 4^ 

4 

A 

Substituting  Equation  (3.28)  into  Equation  (3.24)  yields: 

CD 

P 


G(t)  - 2P  t 
ave 


I1 + S fe)# “ ■ 


(3.28) 


(3.29) 


The  important  feature  of  Equation  (3.29)  is  the  relative  decrease 
in  the  contribution  from  harmonics  as  cot  becomes  large.  It  is  also 
evident  that  the  dominant  contribution  to  temperature  increase  comes 
from  the  average  power  contribution.  Neglecting  the  harmonic  terms 
gives : 

1/2 

G(t)  = 2P  t 9 (3.30) 

ave 

which  is  the  response  to  a square  wave  pulse  of  power  P . If 

1/2 

2P  t * < 2K  (3.31) 

ave  — 


over  many  cycles , then  damage  must  be  assessed  using  the  exponential 
decay  factor.  For  cyclic  responses,  the  power  dissipation  will  be 
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00 


given  by  the  form: 


P(t)  = e (V + P^  sin  nu)tj  = e P(t) 


ave 


n-1 


(3.32) 


where  P(t)  is  given  by  Equation  (3.23)  Based  on  the  analysis  of  the 
undamped  case,  we  assume,  without  rigorous  proof,  that  for  times  long 
enough  so  that  damping  becomes  important,  the  dominant  contribution  to 
G(t)  will  result  from  the  first  term  in  Equation  (3.32),  namely, 


P(t)  = e~Yt  P 


ave 


(3.33) 


Using  the  Equation  (3.33),  the  expression  for  G(t)  becomes: 

t 


G(t) 


ave 


f e~Yt ' 

1 (t-t’)1/2 


dt' 


(3.34) 


Making  the  substitution 


x = Yt  , x'  = Yt ' 


then  yields 


G(t) 


(3.35) 


(3.36) 


Now  letting 


y = x'/x 


(3.37) 


Depending  on  the  relationship  between  the  device  voltage  and  the  input 
voltage,  the  power  damping  rate  for  the  device  may  be  either  y or 
2y.  For  the  purposes  of  this  discussion,  we  use  Y* 
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gives: 


* 


« 


““  / 7$*  “ 

- pavetl/2  2[1  - x / (l-y)1/2  e_MX  dy  ] 

0 

2P  . /0  2P 

= — Jj2  x [l-x<Kx)]  = — y/2  ♦ 

y y 

where 

1 

<f>(x)  - J'  (1-y)1^2  e yX  dy 
0 


<Kx) 


Xl/2  [ 1 - x<J>(x)  ] 


(3.38) 


(3.39) 


(3.40) 


(3.41) 


(3.42) 


Figure  6 shows  a plot  of  <Kx)  vs.  x.  For  small  values  of  x, 
<J>(x)  is  approximated  as: 

1 

<Kx)  * J (l-y)1/2  dy  = | (3.43) 

0 

which  gives  for  x = yt  < 1 : 

2p 

G(t)  " — 172  (yt) 1/2t 1 " f(Yt^  = 2PaVetl/2^ 1 " 3 (Yt)]  <3.44) 

Y 

»-  y X 

On  the  other  hand,  for  large  values  of  x the  e term  in  the 

integrand  will  be  the  controlling  factor  in  the  integration.  In  the 

—TJX 

region  where  e is  non-negligible  y will  be  very  small  compared 
to  unity  and  <f>  will  be  approximately  given  by: 
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i 


4>(x) 


l|>(x) 


Figure  7.  ^(x)  vs.  x 
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« 


(3.45) 


Substituting  Equation  (3.45)  in  Equation  (3.40)  gives  for  Yt  > 1, 


G(t) 


2P 


ave 

1/2 


x1/2[l-(l-e"X)] 


2P  -i  to 

ave  1/  2 -x 

yl/2  X 6 


(3.46) 


Figure  7 shows  a plot  of  the  function  <Kx) , which  is  propor- 
tional to  G.  As  observed,  t|»(x)  reaches  its  maximum  at  the  point 


xsl 


(3.47) 


or  equivalently  at  time 


1 

tss  ~ 

. y 


(3.48) 


The  maximum  value  of  t(>(x)  is  0.55  which  then  gives 


2P 


max 


_ave_  , .1/2 

1/2 


1.5P 


ave 


1/2 


(3.49) 


A 
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4.0  CONCLUSION 

In  this  report,  we  have  developed  several  analytical  techniques 
which  can  be  used  for  evaluating  permanent  damage  of  semiconductor 
devices  which  are  contained  in  interface  circuits.  These  techniques 
are  intended  for  use  in  a hand  analysis  and/or  the  screening  phase  of 
hardness  evaluation. 

Included  in  this  discussion  are:  the  effects  of  finite  size  p-n 

junctions  on  maximum  temperature  buildup;  an  evaluation  of  inter- 
pulse cooling;  and  a determination  of  the  response  of  these  devices 
to  damped  periodic  waveforms.  Limitations  of  the  results  are  presen- 
ted. 
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